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For electron-phonon Hamiltonians with the couplings Hnear in the phonon operators we construct 
a class of unitary transformations that separate the normal modes into two groups. The modes in 
the first group interact with the electronic degrees of freedom directly. The modes in the second 
group interact directly only with the modes in the first group but not with the electronic system. 
We show that for the n-level electronic system the minimum number of modes in the first group is 
Us = (n^ + n — 2)/2. The separation of the normal modes into two groups allows one to develop 
new approximation schemes. We apply one of such schemes to study exitonic relaxation in a model 
semiconducting molecular heterojuction. 



I. INTRODUCTION 

One of the central assumptions in formulating the re- 
laxation dynamics of a quantum system embedded in a 
condensed media is the partitioning between degrees of 
freedom treated explicitly and implicitly. Often, we are 
guided by our "physical intuition" in separating the sys- 
tem (explicit degrees of freedom) from the bath (implicit 
degrees of freedom). This intuitive approach breaks- 
down severely when the coupling between system and 
bath degrees of freedom become too strong. This situ- 
ation is exasperated when one attempts to treat some 
of the degrees of freedom using quantum mechanics and 
the remaining degrees of freedom using classical mechan- 
ics. What is clearly lacking is a unambiguous approach 
for formally partitioning a system with both discrete and 
continuous degrees of freedom into interacting sub-spaces 
such that the motions that are most strongly coupled 
are given explicit treatment while the remainder serve as 
auxiliary degrees of freedom (quantum or classical), as a 
thermal bath, or are ignored completely. Clearly, such 
a partitioning is dependent upon the specific system at 
hand. However, the prescription for making the parti- 
tioning should be universal and depend only upon the 
form of the original Hamiltonian. 

In many respects, we anticipate the presence of "spe- 
cial" degrees of freedom. One of the central tenants of 
modern chemical physics is that of the reaction coordi- 
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nate, espectially in the Marcus-Hush model of charge- 
transfer between donor and acceptor species. The re- 
action coordinate that parameterizes the free-energy 
parabolas in this classic model are not explicitly defined 
in terms of actual nuclear motions. The driving force, 
electronic coupling, and reorganization energy terms ap- 
pearing in the Marcus-Hush expression for the electron- 
transfer rate are all experimentally derived quantities. 
Nuclear motion enters into the description when we make 
the connection between the polarization field and the re- 
organization and relaxation of the donor and acceptor 
molecules and their surroundings. 

In this paper we present a procedure for deducing 
the reaction coordinate for a multidimensional electron 
transfer process. By using a series of canonical trans- 
formations we can uniquely decompose the fully-coupled 
problem into a subset of strongly correlated vibronic mo- 
tions that have clear dynamical and structural meaning 
and carry all of the electron/phonon couplings. This sub- 
space is coupled to the remaining vibrational degrees of 
freedom which in turn are completely decoupled from 
the electron transfer dynamics, but serve as a dissipative 
subsystem. 



II. PROJECTED MODES 



A wide class of the electron-phonon systems is described by the following Hamiltonian 

n n N N 2 ^ 2 2 

H = J2ea\a){a\+ ^ ^g,,,:,„|a)(6| + ^ | + ^ (1) 

a—l a—b—1 i—1 i—1 i—1 

Here |a)'s denote electronic states with energies Ca, Xi and pi are coordinate and momentum operators for the normal 
mode i with frequency Wj, and are the coupling parameters of the electron-phonon interaction. N is the number 
of normal modes and n is the number of electronic states. We seek to rewrite H in terms of the new coordinate Xi 
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and momentum Pi operators so that it has the following form 

a—l a—b—1 i—1 i—lj—ns+1 i—1 i—1 

where rtj — "^'^""'"^ and uj^ are the new oscillator frequencies. In the new Hamiltonian only modes are directly 
coupled to the electronic degrees of freedom with 5^^^ being the new coupling parameters. We will refer to these modes 
as the system modes |3J H] . These Us modes are also coupled to the N — remaining normal modes (that are not 
directly coupled to the electronic degrees of freedom). We will call these N — modes the bath modes. Coefficients 
7j'j give the couplings between the system modes and the bath modes. The transition form Eq. (jlj) to Eq. (j2| can be 
achieved by a suitable linear transformation of the coordinate and momentum operators. 
Consider the following linear transformations 

X, = E^^y-^J' (3) 

j 

P. = J2{M-%Pj. (4) 

Here Mij are the elements of an arbitrary real matrix M. We will assume that M has an inverse, M"^ M"^ denotes 
the transpose of M~^. It can be shown that transformations ([s]) and Q are canonical, i. e., Xi and Pi satisfy the 
usual commutation relations for coordinates and momenta. Since matrix M is real the new operators Xi and Pi are 
Hermitian. The inverse transformations are given by 

i 

= E^^J^J- (6) 

j 

In Appendix A we show that transformations (|3| and Q can always be viewed as unitary transformations of operators 
Xi and Pi. 

Our goal of going to the new coordinates can be achieved with orthogonal matrices M. Therefore, in this section we 
will assume that M is orthogonal, i. e., M = M^^. We will discuss the possible uses of the non-orthogonal matrices 
below. 

In order to obtain the explicit form of M we substitute expressions ^ and ^ for Xi and pi into Eq. (jlJ. This 
gives 

H = Y^e^\a){a\+ E E •^^''^A'^^ l«) (^1 + E ^ + ^ ^^^^^ 

a—l a—b—1 i—1 j — 1 i—1 i—j — l 

Here 

Note that the momentum operators remain uncoupled for any orthogonal M in Eq. ^ because the matrix associated 
with the quadratic form p^ is the unit matrix and is invariant under the orthogonal transformations. 

I 



Hamiltonian Q will have the form ([2| if the following 
requirements are satisfied. First, the summation over j in 
the second term has to extend only from 1 to rig. Second, 
the mode-mode interaction part for the coordinate oper- 
ators must be such that there is no interaction within the 
group of the system modes and no interaction within 
the group of the N — Us bath modes while the interaction 
between these to groups remains. 



Consider the second term of Eq. ([t]) . It contains the 
sum "^ij gabiM~j^Xj. This expression involves multipli- 
cations from the left of the row vector gab with elements 
gabi on the N x N matrix M~^. Since gabi — gbai there 
are ris — n{n — l)/2 linearly independent vectors gab for 
the n level electronic system. To simplify the notation 
we will denote these tt-s independent vectors g^ where a 
can run from 1 to ris- Recall that the columns of any 
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orthogonal matrix form a complete set of orthonormal 
vectors. Then, our first requirement of the reduction of 
the number of operators appearing in the electron- 
phonon interaction to Ug can be viewed as the search for 
orthogonal matrices M^^ whose N — Ug columns are or- 
thogonal to 71s vectors g^. Obviously these matrices are 
not unique. To obtain the columns of such matrices we 
construct Us arbitrary orthonormal vectors in the sub- 
space of rig vectors go. and N — Ug arbitrary orthonormal 
vectors in the subspace that is orthogonal to all go's. We 
will call these two groups and k^, respectively. 

The unique choice of matrix M^^ is achieved by satis- 
fying our second requirement. Explicitly, vectors k^ and 
k^ must be such that the Ug x Ug Hessian matrix for the 
system oscillators 



N 



(9) 



N X N Hessian matrix for the primitive modes following 
the projection 



K = Pcj^P + Qcj^Q. 



(14) 



Both Pw^P and Qw^Q are N x N matrices. However, 
the first will only have Up non-zero eigenvalues {cOg} 
with corresponding eigenvectors {Mg} while the second 
(residual set) will have N — Ug non-zero eigenvalues {uj^} 
with corresponding eigenvectors {Mr}. Thus the full 
N X N transformation matrix is formed by joining the 
non-trivial vectors from the P subspace to the N ~ 2 
vectors from the Q subspace M = {Mg, Mr}- 

The coupling between the V and Q subspaces, j'^a', 
are then given by the non-diagonal blocks of the hessian 



transformed into the eigenbasis of K. 



(15) 



and {N - Ug) x {N 
oscillators 



Ug) Hessian matrix for the bath 



N 



(10) 



are both diagonal. The eigenvalues of these matrices will 
give the squares of the new oscillator frequencies. The 
coupling coefficients 7^^, in Eq. Q are given by 



N 



y' ,^^k' 

laa' / , ""o 



(11) 



The first subscript in 7^^' refers to the system modes 
and the second one to the bath modes. Coefficients g^^^ 
for the electron-phonon coupling in Eq. ([2| are 



9abo 



E 



9abi 



(12) 



An alternative view on the above derivation can be ob- 
tained if we introduce a projection matrix onto the sub- 
space of vectors go. . We do so by defining the projection 
operator 



' g/3 



(13) 



a/3 



with Safi = ga 'g/S and (g) denotes the outer product. The 
summation is restricted to only linearly-independent vec- 
tors. This is an X matrix projects out all primitive 
modes that are directly coupled to the electronic degrees 
of freedom and its complement Q = I — P projects out 
all motions not drectly coupled. It is simple to show that 
P Q = Q P. 

We can now form the matrix M^^ that transforms 
between the primitive and projected coordinates by find- 
ing the eigenvectors of the block-diagonal elements of the 
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(Ms-L-M,), 



(16) 



where Ms is a ris x A^ matrix whose rows are vectors {Mg} 
and Mr is a A^ — ris x A matrix whose rows are vectors 
{Mr} Finally, the electron/phonon coupling coefficients 
9abi obtained as 



Ms 



(17) 



In short, the projection scheme based upon the linear 
coupling coefficients provides a robust and efficient way 
to extract a subset of motions that are strongly coupled 
to the electronic degrees of freedom. 

The final hessian matrix now has the form: 



K 



7 

7 UJr 



(18) 



where Ug are the frequencies of the ng coupling modes, 
LUr are the frequencies of the residual N — rir modes. 
The two types of modes are coupled via 7. One can 
easily verify that this new hessian is related by orthogonal 
transformation to the original (diagonal) hessian matrix 
and obtain two sets of Boson operators Ai G V with 
[As, a);,] = See' and Br e Q with [Br.Bl,] = 5rr' which 
separately act in the V and Q subspaces respectively so 
that [As,i?r'] = 0. Clearly, modes in Q are orthogonal 
but coupled to the modes in V. 

Several generalizations of the transformations consid- 
ered above are possible. Any Hamiltonian in which the 
electron-phonon coupling involves Ug sums of the form 
aiXi or hiPi where a; and hi are constants can be 
transformed in such a way that only Ug modes are di- 
rectly coupled to the electronic degrees of freedom. An 
example of such Hamiltonian can be found in Ref. [T]. 
Note that it is also possible to further reduce the number 
of the system modes by one, if the completeness relation- 
ship for the electronic degrees of freedom is used. So far 
we considered the unitary transformation that isolates 
the minimum number of modes coupled to the electronic 
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system. Imposing different requirements on the unitary 
transformation we can obtain different decompositions of 
the phonons into groups. Certain schemes may be better 
suited to a particular problem than others. For example, 
in Ref. [U [31 Uj the so-called Mori chain transformation 
was applied in the case of two electronic states. Here, 
the primitive Hamiltonian is first transformed as above 
so that only three normal modes are coupled to the two- 
level system and the bath contains A'' — 3 modes. How- 
ever, at the next stage a new transformation is applied 
to the bath modes so that the three system modes are 
coupled directly to only three other bath modes while 
the new bath consists of — 6 modes. The procedure is 
iterated until all the modes are arranged into a hierarchy 
whereby each level contains three modes coupled to the 
three modes in the previous level and three modes in the 
next level. Note that the modes within each level of the 
hierarchy are not coupled direct to each other but only 
indirectly via the neighboring levels of the chain. We 
recently applied this type of transformation to study ex- 
citon relaxation in a conjugated polymers heteroj unction 
Refs. [51|5J[7j and in the photoreactive yellow protein [8J. 
Alternative form of the Mori chain can be constructed. It 
is possible to develop a hierarchy in which modes within 
each level are coupled to each other but each of them is 
coupled to only one mode from the previous level and 
one mode from the next. Without going into details we 
note that to construct such a chain requires using non- 
orthogonal matrices whose rows involve vectors formed 
from coupling coefficients without preliminary orthogo- 
nalization as above. Note, however, that the new oper- 
ators are still related to the old ones through a unitary 
transformation as shown in Appendix A. 

In Appendix B, we show that that it is possible to 
transform the Hamiltonian in Eq. [2] to the form in which 
there are still Ug system modes but each of the operators 
\a){(3\ + 1/3) (a I (for a ^ fi) and |q;)(q;| interacts with its 
own single mode. In particular, for the two-level system 
only one mode is directly responsible for the electronic 
transitions. 



III. ELECTRONIC RELAXATION IN A 
MOLECULAR HETEROJUNCTION DIMER 

As an application of our approach we consider the de- 
cay of an excitonic state into a charge-separated state 
at the interface between two semiconducting polymer 
phases (TFB) and (F8BT). The chemical structures of 
these are given in Fig.[T] Such materials have been exten- 
sively studied for their potential in organic light-emitting 
diodes and organic photovoltaics [9} [T0 | [TT | [T2 l [T3 l [14] . 
At the phase boundary, the material forms a type-II semi- 
conductor heteroj unction with the off-set between the va- 
lence bands of the two materials being only slightly more 
than the binding energy of an exciton placed on either 
the TFB or F8BT polymer. As a result, an exciton on 
the F8BT side will dissociate to form a charge-separated 




FIG. 1: Chemical structures of F8BT and TFB oligomers. 

(exciplex) state at the interface. [HI HH [131 [151 US] 

F8BT* : TFB — > F8BT" : TFB+. 

Ordinarily, such type II systems are best suited for photo- 
voltaic rather than LED applications However, LEDs fab- 
ricated from phase-segregated 50:50 blends of TFB:F8BT 
give remarkably efficient electroluminescence efficiency 
due to secondary exciton formation due the back-reaction 

F8BT" : TFB+ — > F8BT* : TFB. 

as thermal equilbrium between the excitonic and charge- 
transfer states is established. This is evidenced by long- 
time emission, red-shifted relative to the emission from 
the exciplex, accounting for nearly 90% of the integrated 
photo-emission . 

Here we consider only two electronic levels correspond- 
ing to 

\XT) = F8BT* : TFB & \CT) = F8BT" : TFB+. 

We take the vertical energies for these two states as 
ecT = 2.019eV and exT = 2.306 eV. As was shown in 
Ref. [T71 [TSl [HI [50], in poly-fluorene based systems, 
there are essentially two phonon bands that are coupled 
strongly to the electronic degrees of freedom as evidenced 
by their presence as vibronic features in the vibronic 
emission spectra. Within our model these are represented 
by two narrow sine-shaped bands. The first, centered 
about 97.40 cm~^ with a width of 10.16 cm~^, accounts 
for the low-frequency torsional motions while the second, 
higher frequency band centered at 1561.64 cm~^ with 
width 156.30 cm~^ accounts for the C=C ring-stretching 
modes. These primitive modes are entirely localized on 
one chain or the other. The vibronic couplings within 
the model were determined by comparison between the 
Franck-Condon peaks of the predicted and observed spec- 
tra of the system [17]. We studied this model system 
before using several different approaches [H [51 [H] W^ . 
As before, we consider the initial excitonic state as being 
prepared in the \XT) state corresponding to the photoex- 
citation of the F8BT polymer. The full set of parameters 
for the model may be obtained as the supplemental ma- 
terial of Ref. [IJ. 
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FIG. 2: Projection of the three couphng modes onto the prim- 
itive modes for the model heteroj unction system, (a) Projec- 
tion onto low frequency torsional modes and (b) projection on 
to the higher frequency C=C modes included in the model. 
Red, Blue, and Green indicate the contributions from cou- 
pling modes 1,2, and 3 respectively. 



FIG. 3: a.) Energy level diagram showing relevant electronic 
states and corresponding vertical excitation energies, b.) A 
slices through the 3-dimensional Born-Oppenheimer potential 
energy surface for the model heterojunction dimer. Here, is 
fixed at the global minima (corresponding to the CT diabatic 
curve). 



A. Contributions from Coupling Modes 

First, we examine the effect of mode-projection on the 
model Hamiltonian describing the heterojunction. Fol- 
lowing the mode-projection, the 28 primitive modes re- 
duce to three coupling modes with frequencies cOs — 
1587.61cm~^, 1625.07 cm^^, and 1637.22 cm"-'^ respec- 
tively. The projection of these modes onto the primi- 
tive basis is shown in Fig. [2] These projected modes 
involve contributions from both the high and low fre- 
quency phonon bands, although the dominant contri- 
butions come from the higher frequency phonon band. 
Furthermore, since the primitive modes are localized to 
the individual molecular units, we can comment that the 
first coupling mode uji = 1587.61 cm~^ (shown in red in 
Fig. |2| involves contributions from both molecules, while 
the other two modes (shown in blue and green, in Fig. [2] 
respectively) are dominated by local contributions com- 
ing from either the F8BT chain or the TFB chain. 

In Fig. [3] we show some details of the potential energy 
surface generated by our model Hamiltonian. Fig. |3^ 
shows the position of the two relevant electronic states 
corresponding to the emissive excitonic state and the 
weakly emissive exciplex state. Within the model, there 
does exist an intermediate charge-transfer state that is 
energetically close to the XT state. We recently explored 
the and energetics [21] and dynamical |6] implications of 
such dark exciplex states using quantum chemical and 
numerically exact quantum propagation. In the initial 
transfer, the intermediate state (for this system at least) 
plays a minimal role with most of the population being 
transferred directly from the XT to CT. However, it does 
play a role in the back transfer process. In Fig. [3]d we 
give a slice along along one of the projected mode coor- 
dinates {X2) keeping 2:3 fixed at the energy minimum of 
the CT state. The slice where xi is at the CT energy 
minimum indicates that the system is nearly in the "bar- 



rierless" regime for electron transfer as suggested by the 
experimental observations. [S] 

B. Population Relaxation 

We now examine the dynamical ramifications of the 
mode projection scheme for the system following exci- 
tation to the XT state. In Ref. [T] we showed that 
to second-order in the electron/phonon couplings, non- 
Markovian contributions to the population decay can be 
included within a time-convolutionless form of the master 
equation (TCLME) 

^Ml^J2wMt)P,{t)-Y,Wp^{t)P^{t). (19) 

in which the time dependent rates Wai3{t) are given by 

W^„(t) = 2Re f c^T(M„/3A^0„(r))e-^(^^"-^>)^ (20) 
Jo 

where 

{M^pMp^ir)) = Tr (M„^M^„(t)p,,) . (21) 

is the autocorrelation of the off-diagonal coupling opera- 
tor in the Heisenberg representation for the canonical en- 
semble Peg. The explicit form of Capir) = {MapMpa{T)) 
is complicated and lengthy and the reader is referred to 
the original work for its details. Provided that the corre- 
lation function Cap{t) — > at long times, the golden-rule 
rate is given by 

Wf,;- = lim Wp^it). (22) 

At long times when the rates become constant the 
TCLME preserves both positivity and detailed balance. 
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FIG. 4: The real part of the correlation function Ci3(t) for 
the fully-coupled 28 modes model (solid gray) and reduced 3 
(solid black) and 6 (dashes) mode decoupled models. Here we 
use a temperature of T = 298K. 



FIG. 5: The TCLME population of the excited state as a 
function of time for the fully-coupled 28 modes model (solid 
gray) and reduced 3 (solid black) and 6 (dashes) mode decou- 
pled models. T = 298K. 



However, at intermediate times one can obtain nega- 
tive rates depending upon the spectral density of the 
coupling [T]. 

In Fig. g] and Fig. [s] we compare the results of vari- 
ous approximate decoupling scheme with results obtained 
without using the decoupling. As a first approximation, 
we disregard all the bath modes and their couplings to 
the system modes. Thus, only the terms describing elec- 
tronic degrees of freedom and the three c-modes are kept. 
Clearly, this can only be a good approximation if in the 
transformed Hamiltonian (Eq. [2]) the coupling between 
the system modes and bath modes is weak. For the case 
at hand, this proves to be a reasonable approximation 
for a short time dynamics but a poor approximation for 
longer times as the population never relaxes from the ini- 
tial state. This indicates that the coupling between the 
c-modes and residual r modes is sufficiently strong and 
cannot be ignored. 

The picture is completely different if we perform inde- 
pendent unitary transformations within each of the two 
phonon bands of our model. In this case the electronic 
degrees of freedom are coupled to a total of six c-modes, 
three for each phonon band. Within the TCLME formal- 
ism, the computed population decay is virtually indistin- 
guishable from the fully coupled (28-mode) case even out 
to long times. 

In Ref. [3 13 the same model was investigated by 
using different truncations of the Mori chain with the 
time evolution of the multi-dimensional electron/phonon 
wavefunction was propagated using the numerically exact 



multi-configuration time-dependent Hartree (MCTDH) 
method 23, 24 . The results obtained there are qualita- 
tively similar to ours in that too severe of a truncation 
produces no relaxation and that both high and low fre- 
quency contributions are needed to affect the fission of 
the exciton. 



IV. DISCUSSION 

Here we have presented one very straightforward and 
essentially universal procedure for reducing the com- 
plexity of a quantum dynamical simulation of a multi- 
dimensional multi-state system by reducing the number 
of normal modes directly interacting with the electronic 
system. Splitting all the normal modes into two groups 
allows one to apply different levels of theory to each of 
these groups. In particular, one can expect that as far 
as the electronic dynamics goes that the modes that are 
directly coupled to the electronic degrees of freedom are 
more important and showed be treated by more exact 
methods while the residual modes can be treated semi- 
classically as a bath. Such approximations should cer- 
tainly be justified if the system modes are weakly cou- 
pled to the bath modes. We showed that rather then 
using the minimum number of system modes it is possi- 
ble to adjust the number of the system modes to reduce 
their coupling the bath modes. It would be desirable to 
develop an unambiguous procedure for choosing the op- 
timum number of modes that are still few in number and 
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at the same time weakly coupled to the residual modes. 
So far the only unambiguous criteria that we have is that 
doing transformations of the type of Eq. ^ within the 
groups of modes with near-degenerate frequencies gener- 
ally leads to the weak system-bath mode coupling. In 
particular, in case of complete degeneracy, the system 
and bath modes are completely decoupled. 



APPENDIX A: 



Consider the linear transformations 



P, 



(Al) 



(A2) 



where M^- are the elements of an arbitrary real invertible 
matrix M. In this appendix we will show that these 
transformations are equivalent to the following unitary 
transformations of operators xt and pi 



X,, 



Pr 



(A3) 



and, therefore, transformations in Eqs. (A2| are canoni- 
cal even though they arc not orthogonal. 

By virtue of the singular value decomposition theorem 
any invertible real matrix M can be written as 



M = O1DO2, 



(A4) 



where Oi and O2 are orthogonal matrices and D is a pos- 
itive definite diagonal matrix. It is sufficient to show that 
that there is a unitary transformation corresponding to 
each of these matrices. Then the unitary transformation 
corresponding to M will be a product of these transfor- 
mations. Consider the orthogonal matrices first. They 
can have a determinant of either 1 or — 1 . If the determi- 
nant of an orthogonal matrix is positive it corresponds 
rotations in the multi-dimensional coordinate space. Let 
us show that in can be represented by the following uni- 
tary transformation of Xi and pi 



X,, = U-^x,U„ 



P, = U-^p,Uo, 



where 



Uo = exp ( - ^ E Aij{piXj + Xjp.,) 



(A5) 



(A6) 



and A is a real anti-symmetric matrix. The explicit effect 



of the unitary transformation A5 can be calculated by 
introducing an auxiliary operator 

Uo (0 = exp ( - J ^ A,, {p,x, + xjp,)q . ( A7) 



Consider the auxiliary coordinate and momentum oper- 
ators 



Differentiating Xi{£^) and Pi(^) with respect to ^ and us- 
ing the commutations relations between coordinates and 
momenta we obtain 



Solving the last two equations subject to the initial con- 
ditions Xi{0) = Xi and -Pi(O) = pi and bearing in mind 
that X, = Pi = P,{1) we have 



X, 



Since any orthogonal matrix with determinant 1 can be 
written as an exponential of an anti-symmetric matrix 
we can reproduce any orthogonal matrix in Eq. A10| by 
a suitable choice of A. Thus, any orthogonal transforma- 
tion with determinant 1 can be represented by a unitary 
transformation of the form of Eq. ( A5 1. Let us now show 



that there exists a suitable unitary transformation of Xi 
and Pi if the determinant of either of the orthogonal ma- 
trices in Eq. |A4| is negative. Any orthogonal matrix with 
the determinant of —1 can be written as a product of an 
arbitrary orthogonal matrix with the determinant of —1 
and an orthogonal matrix with the determinant of 1 



0_ 



(All) 



The case of the positive determinant orthogonal matrix 
was considered above. Therefore, we can choose a sim- 
ple orthogonal matrix with determinant —1 for 0'_i and 
show that there is a corresponding unitary transforma- 
tion. A simple choice is a diagonal matrix one of whose 
diagonal elements is equal to — 1 and the rest are equal to 
1 . Its action corresponds to the space inversion for one of 
the oscillators and has the following unitary realization 



X, 



U^^X.Uq, 



P, = U-^P,Ug 



with 



Ug — exp iTT 



2 



(A12) 



(A13) 



Finally matrix D being a positive definite diagonal ma- 
trix corresponds to the rescaling of the coordinate and 
momentum operators of the following form 



X,,^e^'x,. 



P, = e-^^p., 



and also has a unitary realization as 

X, = U-^x,Us, P^ = u-'p,u,, 

where 



exp 



■PiXi) 



(A14) 
(A15) 

(A16) 



Thus , we showed that transformations in Eqs. (A2| can 



be viewed as unitary transformations of operators Xi and 
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APPENDIX B: 



3x3 matrix L 



In this appendix we will show that it is possible to 
transform Eq. Q to the form in which there are still 
system modes but each of the operators \a) {b\ + \b) {a\ (for 
b) and \a){a\ interacts with its own single mode. We 
consider the case of a two-level system, generalizations 
to higher-level systems are trivial. Consider the case of 
two-level electronic system and, therefore, three system 
modes (71 <; ~ 3). Let us show that it can be transformed 
into 

2 

H = Y,ea\a){a\+Y,\l){l\+Y,\2){2\ + Y,{\l){2\ + \2){ 

3 N 3 3 



5112 5113 \ 
L = I 5221 5222 5223 I 
\5i2i 5122 5123/ 



(B3) 



and {M~^)ij denote matrix elements of the inverse of the 
transpose of L. Matrix elements Uij and Vij in Eq. Bl 
are given by 



Uij — ^ ^ LkiLkj, 



Y,u',\L-%,{L-%,{m) 



k=l 



k=l 



and 



1=1 j=4 
^ p2 



«j=i 



(B5) 



2 ' ^ 2 ' 

i=4 2=4 



fc=l 



(Bl). 



where we have written the electron-phonon interaction 
part explicitly. Thus, only one mode is directly responsi- 
ble for the electronic transitions. Eq. |Bl| can be obtained 
from Eq. [2] by the following non-orthogonal transforma- 
tion of coordinate and momentum operators in the space 
of system modes 



We would like to stress again here that even though ma- 
trix L is non-orthogonal operators Yi and Ilj are unitar- 
ily equivalent to, respectively, Xi and Pi, and therefore, 
Hamiltonian (Bl) is unitarily equivalent to Hamiltonian 
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This is a special case of the coordinate transformations 
we used earlier. Here Lij are elements of the following 
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